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ABSTRACT 

We consider the evolution of a system containing a population of massive planets 
formed rapidly through a fragmentation process occurring on a scale on the order of 
100 au and a lower mass planet that assembles in a disc on a much longer timescale. 
During the formation phase, the inner planet is kept on a circular orbit due to tidal 
interaction with the disc, while the outer planets undergo dynamical relaxation. In- 
teraction with the massive planets left in the system after the inner planet forms may 
increase the eccentricity of the inner orbit to high values, producing systems similar 
to those observed. 

Key words: giant planet formation - extrasolar planets - dynamical relaxation - 
orbital elements 



1 INTRODUCTION 

The extrasolar planets discovered so far have masses, semi- 
major axes and eccentricities in the range 0.16-11 Jupiter 
masses (Mj), 0.038-4.47 au and 0-0.93, respectively. In gen- 
eral, high eccentricities are very difficult to explain in the 
context of a model where an isolated planet forms in a disc, 
as the disc-protoplanet interaction leads to strong eccentric- 
ity damping (Nelson et al. 2000). 

Papaloizou & Terquem (2001, hereafter PT01) investi- 
gated a scenario in which a population of planetary mass 
objects formed rapidly through a fragmentation process oc- 
curing in a thick disc or protostellar envelope on a scale 
of 100 au. Such a system then underwent dynamical relax- 
ation on a timescale of hundreds of orbits which resulted in 
ejection of most of the objects. It was found that the charac- 
teristics of massive eccentric extrasolar planets and the mas- 
sive 'hot Jupiter' observed so far might be accounted for by 
such a model. However, planets with masses smaller than 
a few Mj are probably too small to have formed through 
a gravitational instability or fragmentation process (Mar 
sunaga & Inutsuka 1999, Boss 2000) and are thus more likely 
to have grown through the core accretion mechanism in a 
protoplanetary disc (Mizuno 1980). 

Here, we investigate a scenario in which we have a 
planet accumulated in a disc (PAID) on a timescale of 
10 6 years, being the observed lifetime of protostellar discs 
(Haisch, Lada & Lada 2001), together with a population of 



massive outer planets already formed through fragmentation 
processes on a much shorter timescale. During its formation 
phase, the PAID is kept in a circular orbit by tidal interac- 
tion with the disc while the system of outer planets under- 
goes dynamical relaxation. After disc dispersal occurring at 
t = 10 6 years, the eccentricity of the PAID can be pumped 
up to high values by interaction with the remaining bound 
massive planets. 

In order to focus on the interaction mechanisms, we 
analyse the motion of a planet under a distant perturber in 
a highly eccentric orbit in § ^[ In § [] we present numerical 
simulations of the evolution of a system containing outer 
massive planets and an inner PAID interacting with each 
other. Finally, in § H we discuss our results. 



2 SECULAR PERTURBATION THEORY AND 
LONG TERM CYCLES 

We consider the evolution of the orbit of a PAID around 
a central star due to gravitational interaction with a mas- 
sive long period perturber on a highly eccentric orbit that 
may have resulted from a prior relaxation process. We il- 
lustrate three possible effects. The first is the appearance 
of long period cycles in which the inner planet may attain 
high eccentricity. The second, occuring for stronger pertur- 
bations, is a sequence of dynamical interactions leading to 
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the collision of the inner planet with the central star. The 
third, occurring for weakly bound perturbers, is an interac- 
tion leading to the gradual ejection of the the outer planet 
to large radii, where additional perturbations from external 
objects could cause it to become unbound, leaving the inner 
planet on a highly eccentric orbit. 

To discuss these three scenarios, we denote the osculat- 
ing semi-major axes of the inner planet and the perturber 
by a and a p , respectively. The corresponding eccentricities 
are e and e p and the masses are m and m p , respectively. The 
central mass is M*. 

When the pericentre distance of the perturber, d p = 
a p (l — e p ), significantly exceeds a, the large difference in 
the orbital periods reduces the significance of short term 
variations such that the interaction can be discussed using 
secular perturbation theory, in which one considers evolution 
of the time averaged orbits. 

Consider first the case when the orbits are coplanar. 
Then analytic treatment is possible in the limit m <C m p . We 
denote (r, cj>) and (r p ,(j)p) the cylindrical coordinates of the 
inner and outer mass, respectively, in a frame centered on 
the star. The perturbing potential energy or Hamiltonian of 
the inner planet, including the indirect term (which accounts 
for the acceleration of the origin of the coordinate system), 
is: 
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For r/r p < 1, we expand H in spher- 



ical harmonics retaining terms up to third order in r/r 
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where P n is the Legendre polynomial of order n. We have 
omitted the term —Gmm p /r p as it does not depend on the 
coordinates of the inner planet (H is defined within the addi- 
tion of a constant). On performing time averages over both 
orbits (see, e.g., Roy 1978), one obtains correct to second 
order in e but with no restriction on e„: 
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Here zu is the angle between the apsidal lines of the two or- 
bits. For m <C m p , we can consider that the orbit of the outer 
planet is fixed, i.e. that a p and e p are constant. The secular 
evolution is characterized by a constant a and oscillations of 
e (see, e.g., Murray & Dermott 1999). Since de/dt tx dH/dzu 
(Lagrange's equation), e is maximum for zu = 0. During 
the evolution H(e,zu) is constant. Since e passes through 
zero (the inner orbit has no initial eccentricity), the max- 
imum value of e, e ma x, can then be calculated by writing 
H(0, zu) = iJ(e ma x,0). This gives e max = 5ae p /[2a p (l - e^)]. 
For illustrative example we take e p = 0.9 and a/a p = 0.01, 
so that d p I a — 10. Then e m ax = 0.118 which is significant 
and independent of m p . 

Below we consider the case a = 0.5 au, a p — 50 au, 
m = 0.3 Mj and m p = 7 Mj. It turns out that for such 
small a, relativistic effects have to be taken into account. 
This can be done following the procedure given in Lin et al. 
(2000), according to which the central potential is modified 
such that: 
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Figure 1. Eccentricities versus time for a/a p = 0.01, m p = 7 Mj 
and m = 0.3 Mj, in the coplanar case (inner planet, short-dahed 
line; outer planet, solid line) and in the non coplanar case (inner 
planet, dotted line; outer planet, solid line), where there is much 
larger variations. In both cases the eccentricity of the outermost 
massive planet remains near to 0.9. 
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This amounts to adding the term -3(GM, 
e 2 /2)/(a 2 c 2 ) on the right hand side of equation 



2 m(l + 
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this one readily finds that e ma x is reduced by a factor 
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parameters and M 4 = 1 Mq, this factor is 1.93, giving 
e ma x = 0.06. 

We have calculated the orbital evolution numerically us- 
ing the Bullirsh-Stoer method (e.g., Press et al. 1993). The 
outer planet was started at pericentre and the inner one on a 
circular orbit at random phase. Of course, in contrast to the 
above analysis, both orbits were allowed to vary. Figure [I] 
shows the evolution of the orbital eccentricities as a function 
of time. The results are in good agreement with the analysis, 
exhibiting eccentricity cyles with amplitide e m ax ~ 0.06 as 
expected. This is much lower than typical values observed 
for extrasolar planets. However, much higher eccentricities 
can be obtained through the Kozai mechanism if the orbits 
are mutually inclined (e.g. Lin et al. 2000 and references 
therein). Analysis for an outer circular orbit suggest high 
eccentricities ma y be obtained for mutual inclinations ex- 
ceeding arccos y/3/5 ~ 40° . We have also followed the or- 
bital evolution of the system when the orbits are initially 
mutually inclined at an angle of 60° and the results are also 
shown in figure |l|. In this case, large amplitude cyles with 
e attaining values up to 0.7 are obtained. This is typical 
for runs of this kind. Cases with a/a p — 0.01 and e p = 0.9 
lead to long term secular variations that can provide high 
eccentricities of the inner orbit. 

When a/d p is larger, the ensuing interaction with the 
inner planet is stronger and can lead to very high eccentrici- 
ties. In such cases, there can be a close approach or collision 
with the central star. An example of a run of this kind is 
illustrated in figure |2| In this case the initial mutual orbital 
inclination, the initial values of a and e p and planet masses 
were as above, but the initial value of a p was reduced by a 
factor of four. The outer planet was started at apocentre and 
the inner one on a circular orbit at random phase. This case 
led to very strong interactions at pericentre which resulted 
in the eccentricity of the inner planet approaching unity and 
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m = 0.3 Mj and initual mutual orbital inclination of 60°. The 
eccentricity of the inner planet {dotted line) approaches unity 
after about t = 7 X 10 4 y. Subsequently, a collision with the 
central star occurs. The eccentricity of the outermost massive 
planet (solid line) remains near to 0.9. 
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Figure 3. Eccentricities (inner planet, dotted line; outer planet, 
long-dashed line) and semi-major axes (inner planet, short- 
dashed line; outer planet, solid line) versus time in the case of 
a weakly bound outer planet. Here a/a v = 0.004, m p = 7 Mj, 
m = 0.3 Mj and the initual mutual orbital inclination is 60°. 
In this case, a decreases slightly as energy transfers to the outer 
planet cause it to reach radii exceeding 10 4 au. At this stage, e 
reaches ~ 0.7, whereas e p is always near to unity. 



a collision with the central star after about t — 7 x 10 y. 
This situation corresponds to a more strongly bound outer 
planet than previously. 

When more weakly bound outer planets are considered, 
the interactions can also lead to high eccentricities of the in- 
ner planet. As an example we consider a case with e p = 0.99 
and a p — 252a, other parameters being as above. The evo- 
lution is illustrated in figure [| Here, d p /a = 2.52, which is 
similar than in the previous case. However, here, the outer 
planet is very weakly bound. Therefore, as a result of sig- 
nificant energy changes occurring at pericentre passage in a 
chaotic manner, it eventually reaches radii exceeding 10 4 au, 
where it could become unbound if additional perturbations 
from external objects occured. The inner planet is left on a 
highly eccentric orbit with e ~ 0.7. 



3 N PLANETS SIMULATIONS 

We consider a system of outer planets each having a mass 
m p — 7-8 Mj, one inner planet with a mass m, and a pri- 
mary star with a mass M» = 1 Mq, moving under their grav- 
itational attraction. We take m = 0.3 Mj (Saturn mass) as, 
among the observed systems, this is the most extreme case 
of a planet forming in a disc. We suppose the outer mas- 
sive planets have formed rapidly through fragmentation on 
a scale of 100 au, while the PAID forms in the disc surround- 
ing the star on a much longer timescale. 



3.1 Initial conditions 

At t = 0, we place 5 < N < 20 massive planets taken at 
random from a mass distribution corresponding to a uni- 
form density spherical shell with outer radius of 100 au and 
inner radius of 10 au. The planets, the orbits of which are 
not coplanar, are given the local circular velocity in the az- 
imuthal direction. As shown by PT01, adopting different 
initial conditions, such as a mass distribution in the form 
of a flattened disc, leads to the same qualitative evolution 
as long as three dimensional relaxation occurs , so our re- 
sults do not depend on these initial conditions. The system 
is then allowed to evolve under the influence of the gravita- 
tional attraction of the central star and the different planets 
(thepotential of the star has been modified according to 
eq. J4|). We also include the tidal interaction between the 
star (which radius is taken to be 1 Rq) and a closely ap- 
proaching planet (see PT01 for details). The approximation 
used is valid only when the planet approaches the star on a 
highly eccentric orbit, which is the case when tidal interac- 
tions first occur in the systems we consider here. The system 
undergoes dynamical relaxation as described by PT01. 

We assume a PAID forms on a timescale of 10 6 y. Dur- 
ing the formation process, any eccentricity that would be 
excited by the outer massive planets is damped by the tidal 
interaction with the disc. We therefore built up the planet 
progressively and constrained its orbit to remain circular 
during this phase, with an orbital radius a in the range 0.3- 
10 au. Depending on a, the planet may have either formed 
at this location or migrated there as a result of tidal inter- 
actions with the disc. 



3.2 Results 

We have run about 30 cases with the initial conditions de- 
scribed above. As found in PT01, most of the dynamical 
relaxation of the system of outer planets is over after a few 
hundred orbits, i.e. a few 10 5 y. This relaxation results in 
most of the planets being ejected, while a few (usually be- 
tween and 3) of them end up on highly, sometimes close, 
eccentric and inclined orbits around the star. These objects 
may still interact with each other after most of the relax- 
ation has occured. As a result, or directly because of the 
main relaxation, they may end up close enough to the star 
to perturb significantly the PAID after it forms. 

When there are still a few outer massive planets left at 
t = 10 6 y, if they make frequent incursions into the inner 
parts of the system, the eccentricity of the PAID is pumped 
up to high values and the object eventually hits the star 
(e.g. the calculation illustrated in figure ^). 
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Figure 4. Semi-major axes (in au, upper plot) and periccntrc 
distances (in au, lower plot) of the N = 4 outer planets (solid 
lines) and the PAID (dotted lines) versus time (in years). Here 
m p = 8 Mj and a = 0.2 au. In this and other similar figures, a line 
terminates just prior to the escape of a planet. Due to interactions 
with the outer planets, the orbit of the PAID acquires a high 
eccentricity, and subsequently begins to be circularized by the 
central star. During this process, the semi— major axis decreases 
whereas the pericentre distance stays constant equal to 0.02 au. 



In some cases, the PAID may approach the star closely 
enough for tidal effects to become important while a colli- 
sion is avoided. This is illustrated in figure ^, which shows 
the evolution of the semi-major axes and the pericentre 
distances of the planets. Here m p = 8 Mj, iV = 4 and 
a = 0.2 au. The interaction with the outer planets left in the 
system at t = 10 6 y increases the eccentricity of the PAID 
up to high values. It approaches the star closely enough for 
tidal circularization to start to occur. This causes the semi- 
major axis to decrease while maintaining the pericentre dis- 
tance almost constant equal to 0.02 au. At this stage the 
inclination between the orbits of the two planets is close to 
90°. 

If they are at most a few incursions down to small radii 
of the massive planets left in the system at t = 10 6 y, the 
PAID may remain on an orbit with at most a moderate 
eccentricity and a constant semi-major axis more or less 
equal to its initial value. This case is illustrated in figure |E| 
which shows the semi-major axes and the eccentricities of 
the different planets as a function of time. Here m p = 8 Mj, 
N = 9 and a = 0.3 au. At t = 10 6 years, there are still 
two massive outer planets with high semi-major axes and 
eccentricities interacting with each other and one massive 
planet on an eccentric orbit with a semi-major axis of about 
10 au. Close approaches of the two outermost planets (which 
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Figure 5. Semi-major axes (in au, upper plot) and eccentricities 
(lower plot) of the N = 9 outer planets (solid lines) and the 
PAID (dotted lines) versus time (in years). Here m p = 8 Mj and 
a = 0.3 au. The PAID acquires an eccentricity which varies in 
a cycle between 0.1 and 0.24 whereas its semi-major axis stays 
constant equal to the initial value. 



eventually are ejected) pump the eccentricity of the PAID 
up. This eccentricity subsequently goes into a cycle where 
it varies between 0.1 and 0.24 due to the secular interaction 
with the third outer object. This is the same type of process 
as illustrated in figure The semi-major axis of the PAID 
stays constant equal to its intial value a = 0.3 au and the 
inclination between the orbits of the two planets oscillates 
between and 30°. 

When the few outer planets left from the relaxation pro- 
cess make at most a few incursions in the inner system, the 
eccentricity of the PAID may be increased up to high values 
by the secular interaction with one of the massive objects 
with appropriate orbital parameters, but without the pro- 
cess resulting in a collision with the star. The eccentricity 
then varies in a cycle as described in § ^. This also happens 
when there is only one outer planet left at t = 10 6 y, provided 
it is close enough to perturb the PAID. This case is illus- 
trated in figure pi Here m p = 7 Mj , N = 29 and a — 0.77 au. 
Among the massive planets, we have represented the only 
one which was not ejected after about 3 x 10 s y. We note 
that the pericentre distance of the outer planet is about 
9 au. These results are then in agreement with those of § ^J, 
which showed similar behaviour of the eccentricity for such 
a ratio of outer pericentre distance to inner semi-major axis 
when the initial inclination was 60° ( see figure]]]). Here the 
inclination between the orbits of the two planets oscillates 
between and 60°. 

A complete survey of all the runs we have performed in- 
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Figure 6. Same as figure ^, but for N = 29, m p = 7 Mj and 
a = 0.77 au. Among the massive planets, only the one which 
was not ejected after about 3 X 10 5 years is represented. Due to 
the secular interaction with the massive planet, the orbit of the 
PAID acquires an eccentricity which varies in a cycle between 
0.05 and 0.88, whereas its semi— major axis stays constant equal 
to the initial value. 

dicates that due to the randomization of the inclinations of 
the orbits by dynamical relaxation in three dimensions high 
inclinations always occur and an eccentric PAID is produced 
provided the semi-major axis of its orbit is in the appropri- 
ate range, that is not less than about 10 times smaller than 
the pericentre distance of the innermost outer giant planet 
and there need be only one of these. If this condition is sat- 
isfied, efficient eccentricity pumping occurs because of the 
inclination of the orbit of the PAID with respect to that of 
the outer planet. An assessment of the probability for the 
semi-major axis of the PAID to be in the appropriate range 
depends on as yet poorly understood mechanisms of planet 
formation and migration in the disc in the presence of per- 
turbation by an outer giant planet. More information about 
this will be obtained with increasing planet detections. We 
have assumed that the disc is removed after 10 6 y. However, 
similar results are obtained if this occurs at a later time 
as long as the PAID orbit is circularized when the disc is 
present. This is the case as long as the disc is removed after 
the major relaxation (see also § |^) . 



4 SUMMARY AND DISCUSSION 

The analysis and results presented above show that the ec- 
centricity of an inner PAID may be increased up to high val- 
ues by secular interaction with a more distant massive planet 
on an eccentric orbit. The mechanism is efficient when the 



ratio of outer pericentre distance to inner semi-major axis 
is about 10 and the orbits of the planets are significantly 
inclined. 

Such a situation arises when the dynamical relaxation of 
a system of outer massive planets takes place on a timescale 
shorter than the timescale for assembling the PAID in a disc. 
After disc dispersal, as was shown by PT01, the PAID would 
be expected to mostly interact with one more massive object 
orbiting the star with semi-major axis on the order of tens 
of au. If the relaxation is not over by the time the PAID 
forms, close approaches with the massive outer planets tend 
to result in the object colliding with the star (see figure ^) . 

Since giant planet formation occurs on a timescale not 
longer than the disc lifetime of ~ 10 6 y (Haisch et al. 2001), 
assuming formation of massive planets occurs on a scale of 
100 au, an eccentric PAID at less than 1 au from the central 
star is more likely to be produced when the number of outer 
planets in the system is N ~ 15, although N = 5 may also 
produce such a system. 

To illustrate the operation of the processes described 
here, we consider observed planets with masses in the Saturn 
mass range. Two of them (HD 16141 and HD 6434) have or- 
bital eccentricities of 0.28 and 0.3, respectively. These mod- 
est eccentricity systems could be similar to that represented 
in figure |. Two others (HD 83443 and HD 108147) have or- 
bital eccentricities of 0.42 and 0.58 respectively. These larger 
eccentricity systems systems could be similar to that shown 
in figure M. In these cases the eccentricities are driven by 
secular interactions with a more closely approaching outer 
planet. Note however that the probability of getting a very 
high eccentricity on a close orbit is reduced because the 
chance of hitting the central star increases as the semi- 
major axis decreases. Nonetheless, in some cases, the in- 
teractions (secular or not) with the outer planets increase 
the eccentricity of the PAID up to very high values without 
the object colliding with the star. The orbit can then be 
tidally circularized by the central star, which leads to a 'hot 
PAID'. Six 'hot Saturns' have been observed so far, with 
0.24Mj < msini < 0.52Afj and 0.038au < a < 0.066au. 
Although such a process is expected to be rare within the 
context of this model, observational bias against observing 
low mass planets at larger distances makes it difficult to as- 
sess to what extent an an alternative mechanism such as 
direct migration is required. 

If the scenario we have studied in this Letter has opper- 
ated, there should be a massive planet with a on the order of 
tens of au and a high eccentricity associated with the eccen- 
tric PAID on a close orbit. In the case represented in figure ^ 
for instance, the amplitude of the stellar radial velocity in- 
duced by the outer planet is about 60 m s" 1 (that of the in- 
ner planet is between 10 and 20 ms" 1 ). Over a short period 
of time (~ 1 y), only the trajectory of the inner planet would 
be resolved, with the mean value of the stellar radial velocity 
drifting by an amount which depends on whether the outer 
companion is near apocentre or not, and which varies be- 
tween a few m s _1 and tens of m s . Twelve systems have 
been selected by Fischer et al. (2001) as candidates for which 
such residual velocity drifts are potentially observable, and 
their results indicate that the kind of systems we have been 
describing here are on the verge of detectability. Among the 
systems selected by Fischer et al. (2001) with a residual drift 
of 50 m s" 1 over a two year period, HD 38529 has a planet 
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with msini = 0.76 Mj, a — 0.13 au and e = 0.27. This is 
similar to the system illustrated in figure H. 
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